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Correlations between particles can lead to subtle and sometimes counterintuitive phenomena. We 
analyze one such case, occurring during the sudden expansion of fermions in a lattice when the 
initial state has a strong admixture of double occupancies. We promote the notion of quantum 
distillation: during the expansion, and in the case of strongly repulsive interactions, doublons group 
together, forming a nearly ideal band insulator, which is metastable with a low entropy. We propose 
that this effect could be used for coohng purposes in experiments with two-component Fermi gases. 



One of the most exciting features about ultracold 
atom gases is the possibility of experimentally studying 
strongly correlated systems with time-dependent interac- 
tions and in out-of-equilibrium situations (for a review, 
see [1]). For instance, a time-dependent tuning of param- 
eters has been used to illustrate the collapse and revival 
of coherence properties of bosons in an optical lattice sug- 
gested to be well described by the Bose-Hubbard model 
[2]. Also, the question of the relaxation of strongly cor- 
related systems initially prepared in a high-energy state 
that typically is not an eigenstate has been studied [3] , as 
well as the intimately related question of thermalization 
[4] . A third context is the expansion of originally trapped 
particles into an empty optical lattice or a waveguide. 

While expansions after turning off all trapping poten- 
tials and lattices are commonly used to measure the mo- 
mentum distribution of the originally trapped system [1] , 
recently, several studies have investigated the expansion 
in one dimensional (ID) geometries with and without 
the presence of an optical lattice along the ID tubes. 
In ID cases, in contrast to the usual time-of-flight ex- 
periments, interactions play a fundamental role. The 
findings include phenomena such as the fermionization of 
the momentum distribution function of hard-core bosons 
and anyons [5] , the asymptotic transformation of generic 
Lieb-Liniger wave- functions to a Tonks- Girardeau struc- 
ture [6] , the emergence of quasi-coherence in bosonic [7] 
and fermionic systems [8] , or disorder- induced effects [9] . 
Such expansions have already experimentally been real- 
ized in ID ultracold atomic gases [10]. 

While strongly correlated phases of bosons in opti- 
cal lattices have been observed and studied in many ex- 
periments, reaching the quantum degenerate regime of 
fermions is more difficult [1, 11]. Only recently, evidence 
for the observation of the Mott-insulator (MI) transition 
of a two-component Fermi gas in optical lattices has been 
presented [12]. Unfortunately, the temperatures in those 



experiments are still too high to observe the expected 
low-temperature insulating antiferromagnetic phase, and 
to make contact with the regime of interest to high-Tc 
superconductivity [13]. Hence, intense efforts are cur- 
rently under way to develop efficient cooling schemes for 
fermions in optical lattices [14]. 

In this Rapid Communication, we study the expansion 
of a two-component Fermi gas in an optical lattice from 
an initial state with a strong admixture of doubly occu- 
pied sites. Our main result is the observation that double 
occupancies (doublons) are dynamically separated from 
the rest of the system and that they group together into a 
metastable state, which is very close to a Fock state. The 
condition for this to happen is that interaction energies 
need to be much larger than the kinetic energy. Our re- 
sults are in qualitative agreement with recent studies that 
have argued, invoking a similar reasoning, that relaxation 
times can be rather long, and one may thus encounter 
metastable states (see, e.g., Ref. [15]). Experimentally, 
this phenomenon has been observed as repulsively bound 
pairs in the case of bosons [16]. We will also discuss the 
potential of our findings as a cooling technique where one 
dynamically generates a low-entropy region in a system 
with arbitrarily large values of the on-site repulsion. 

Model and methods - We consider the one-dimensional 
(ID) Hubbard model: 

L-l L 

Ho = -tY^ {cl+i^^c,^^ + h.e.) + UY^ n,^in,,i . (1) 

i=l 1=1 

Standard definitions are employed in Eq. (1) [8]. Open 
boundary conditions are imposed, L is the number of 
sites and N the number of particles (with a filling fac- 
tor n = N/L). The nonequilibrium dynamics is stud- 
ied using the adaptive time-dependent density matrix 
renormalization group method (tDMRG) [17] and time- 
dependent exact diagonalization (ED) techniques [18]. 
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FIG. 1: (Color online) Expansion from an initial state with 
Jiinit = 1.8 (iV = 18, U = my. (a) Density {ni(r)) and 
(c) local charge fluctuations Srii for times r f = 0, 10, 20, 30. 
Insets: (b): (ni(r)) vs time for i = 1, 2, 3. (d) Average double 
occupancy d^(T)/x for x = 3,5,10 (dxir) = J2i=o(d'i('^)))- 
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FIG. 2: (Color online) Average double occupancy d3(r)/3 
(?^init = 1.8): vs (a) time r and (b) vs rt'^/U for U/t = 
4, 8, 20, 40, 100 (tDMRG). Insets in (b): di for iV = 5 {U/t = 
40, 100); d2/2 for Af = 3 ((7 = AOt, thin dashed line) (ED). 



We prepare initial states with particles confined into a 
finite region of an optical lattice, with a filhng factor of 
1 < ninit < 2 in that region and zero otherwise. To this 
end [8] , we apply onsite energies to only a portion of the 

system: iJconf = Si^i i^i ~ lO*"^ for 1 < zq < i < £ 
and ei = elsewhere). Hence, at time r < 0, we have 
H = Hq + Hconi, while we turn off i?conf at r = 0'*'. In 
our tDMRG runs, we use either a third-order Trotter- 
Suzuki time-evolution scheme or a Krylov-space based 
method [18], with time steps of Srt = 0.01 and 0.005 
(we set h = 1). During the time-evolution, we keep up 
to 1600 states. 

Results - Figure 1(a) shows snapshots of density pro- 
files {ni{T)) at several times for U = AQt and riinit = 1.8 
(for results for other njnit, see [19]). Since U ^ 4t, single- 
particle hopping can only propagate a particle out of a 
doubly occupied site into the empty sites to the right of 
the initially occupied region if, at the same time, the re- 
maining particle is moved to the left into a site that was 
previously singly occupied, in order to preserve the to- 
tal energy. For that reason, the density in the first sites 
increases as the remaining of the block slowly "melts" . 
From Fig. 1(b), one can see that the density in the first 
sites gets very close to {rii) = 2, which promotes the pic- 
ture that dynamically, doublons are separated from the 
rest and group together in a region of the lattice. 

We introduce the term quantum distillation for this 
process. Figure 1(c) shows the local charge fluctuations 
Snf = (nf) — (rii)'^. Consistent with the picture described 
before, these fluctuations arc the largest in the interface 
region (sites 7-11 in the figure), while they die out in 
the first sites as time increases. To gain a better un- 
derstanding of how this happens with time, we compute 
the average double occupancy d^/x {dx = X]i=o('^« ('''))' 
di = riifnij) on the x first sites, counting from the left. 
Figure 1(d) depicts this quantity for different values of x. 
For a; < 10 and, e.g., a; = 3,5, the average double occu- 
pancy increases towards d^/x = 1, which is quite stable 



over a long time window for U > AOt [Fig. 1(d)]. For 
X = 10, i.e., the block with (n;) > at r = 0, this quan- 
tity remains almost constant: this is basically equivalent 
to saying that only a small portion of the interaction en- 
ergy i?int = converted into kinetic energy 
during the times simulated. 

Time scales - We next address the question on which 
time scales (i) the Fock state is formed and (ii) it decays. 
The latter is equivalent to asking on what time scale a 
band insulator (BI) in the large U limit would decay. 
While from energy considerations it is clear that a piece 
of doubly occupied sites, i.e., a BI, delocalizes on a time- 
scale of U/t'^ [16], we are primarily interested in the case 
of 1 < riinit < 2. To address this situation, we first 
perform ED calculations for a short initial state with 1 < 
"iiiit < 2, specifically iV = 3,5 on two and three sites, 
respectively [see the insets of Fig. 2(b)]. In this case, one 
can clearly see that doublons form on the first sites on a 
?7-independent time-scale, governed by the bare hopping 
matrix element t. Doublons then get dclocalizcd over the 
neighboring empty sites on a time-scale U /t^ (times a 
prefactor that scales with N), consistent with the decay- 
rate of a BI. 

Our tDMRG results presented in Figs. 2(a) and (b) 
unveil that the same picture holds for n = 1.8 with N = 
18. Figure 2(a) shows the average double occupancy on 
the first three sites for several values of U . Quantum 
distillation is best realized for U > 8t (as opposed to 
U < 8t), since, at [/ = 8t, the average double occupancy 
increases with time (r > 5/t), but it reaches a maximum 
value that is well below one. Nevertheless, as the largc-C/ 
results show, the formation of a quasi- Fock state happens 
over a J7- independent time-scale. Figure 2(b) displays 
the same data versus tP/U, and the small-[/ curves fall 
on top of each other at times r t^/J7 > 4, while the large- 
U data seem to approach the same curve, too. The results 
of Fig. 2(b) show, in the example of x = 3, that the Fock 
state will start to delocalize after a time r ~ 50/t in the 
U = 40t case (r ~ 120/t for U = lOOt). 
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FIG. 3: (Color online) Entanglement entropy Sa; (ninit = 1.8, 
U = 40t). Inset: for = 3, 10 at [/ = 4t. 



While wc have presented results for N = 18, the 
physics does not qualitatively change upon increasing N 
at a fixed Tiinit [19]. Essentially, the time scale for build- 
ing up the Fock state scales linearly with N [19]. Note 
that in typical experiments with ID optical lattices, on 
the order of 100 atoms are confined in one ID structure 
[1]. It should be readily possible to experimentally ob- 
serve the quantum distillation by measuring, e.g., the ra- 
dius Rd of the double occupancies [i?^(T) oc '^i{di{T))i'^ , 
see [19]]. We predict that, for C/ ^ 4t and ninit > 1, this 
quantity should decrease as a function of time. As U jt 
depends exponentially on the lattice depth Vb, one can 
enter into the regime C/ ^ 4t by tuning Vq [20] . 

Entanglement entropy - Further quantitative informa- 
tion can be gained by invoking the entanglement entropy 
Sx = —tr[p\n{p)] (Ref. [21]), where p is the reduced den- 
sity matrix of a subsystem of length x [22]. The sub- 
systems of primary interest here are those where dou- 
ble occupancy increases as time evolves, as illustrated in 
Fig. 1(d) (x counts the sites starting from the left end). 

Figure 3 depicts one of our most important results, 
constituting a defining property of the quantum distilla- 
tion process, namely, the spontaneous reduction of the 
entanglement entropy in the metastablc region where 
doublons group together. We plot Sx for several x < 10 
at t/ = AOt, which ought to be contrasted against the 
corresponding U = 4t data, shown in the inset. Fig- 
ure 3 shows that for the leftmost sites, Sx{t = 0) is 
nonzero and remains approximately constant for some x- 
dependent time window {e.g., rt < 5 for a; = 3). At later 
times, the behavior strongly depends on the value oiU/t. 
For U/t = A, S3 increases as the density decreases, while 
for U/t = AO, S3 decreases, on the time scales simulated, 
by a factor of 15 as the metastable state with (ui) =2 is 
generated. We relate this decrease of Sx to the quantum 
distillation, which creates regions with low entanglement. 
This is complementary to experiments aiming at prepar- 
ing maximally entangled states [23]. 

Harmonic trap - While there are experimental efforts 
directed at engineering box-like traps (see, e.g., Ref. [24]), 
which is the case our results presented so far directly ap- 
ply to, in most typical experiments with optical lattices, 
the atoms experience a harmonic confinement [1]. We 



now consider this case to show that quantum distillation 
works even in the more difficult case of strongly inhomo- 
geneous initial states, using -f/^conf — ^rap 1(^ — 
where Vtrap is the curvature of the harmonic potential. 
The confinement gives rise to a shell structure in the 
fermionic density [25]: metallic and MI shells alternate, 
depending on the characteristic density p = N y^Vtiap/t 
[25]. In order to achieve optimal distillation and thus the 
formation of a low-entropy region, one does not need to 
follow the approach often discussed in the literature (see, 
e.g., Ref. [14]), in which U/t is increased adiabatically. 
Instead, we propose that one can start with a trapped 
system with a low value oi U/t = 2. In such cases, dou- 
ble occupancy is energetically favorable against trapping 
energy and hence ideal for our distillation scheme, and 
moreover, the presence of MI shells is suppressed. At 
time T = 0, one can then quench U/t to a large value 
[20] and turn off the trapping potential, but keeping the 
optical lattice on. In many experiments, the trapping 
potential is provided by the same lasers that produce the 
lattice, which, however, is not necessary. In order for our 
scheme to work, one needs a trapping potential that can 
be controlled independently of the lattice [10]. 

The tDMRG results for the time-evolution of the 
quench described above are displayed in Fig. 4. We 
have taken the final values oi U/t after the quench to 
heU = m [Fig. 4(a)] and U = Unit [Fig. 4(b)]. Both 
figures show the average double occupancy dx{T)/x for 
several values of x. In Fig. 4(a), one can see that dx{T)/x 
first decreases and then increases to form the quasi-Fock 
state. It is also illustrative to compare the initial density 
profile to one at later times [inset of Fig. 4(b)]: while 
fast particles escape, a block of particles remains in the 
initial region of the system. A very different behavior of 
dx can be seen in Fig. 4(b). There, since the final value 
of U /t is small, the system simply melts and no signature 
of quantum distillation is seen. 

During the expansion, the entanglement growth [26] 
induced by the quench from U = 2< to AOt competes 
with the reduction of the entanglement entropy Sx of 
blocks with an increasing double occupancy due to the 
quantum distillation. It is thus an amazing result that, at 
sufficiently long times and despite the large final value of 
U chosen in the quench, the quantum distillation wins, as 
is shown in Fig. 4(c). There, we display Sx for x = 2, 3, 4, 
and indeed, these Sx first increase due to the quench, but 
eventually drop below their initial value (dotted lines) 
at the maximum times simulated. For the time-scales 
considered in Fig. 4(c), wc find that S2{t) < 0.25S'2(t = 
0). The same behavior emerges for other values of p, 
Unit, and U [19]. 

So far, our results suggest that quantum distillation 
is robust against a variation of initial conditions and 
quenches. We envision it as a way to experimentally 
achieve very low temperatures in fermionic systems by 
creating a very low entropy band insulating state with 
arbitrarily large values of U/t. Once such a state is cre- 
ated, it can be used as an initial state to achieve an an- 
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that trapping potential and of the ratio U/t. That way 
one could produce a final low temperature Mott insulat- 
ing state starting from a BI. 

In conclusion, we demonstrated that low-entropy 
states of two-component Fermi gases can be dynamically 
created in optical lattices utilizing the expansion of par- 
ticles into the empty lattice in the limit of strong interac- 
tions. After having the low entropy state in an appropri- 
ate trap, an adiabatic lowering of trapping potential and 
interaction strength can lead to a final low temperature 
MI. We stress that the quantum distillation we discussed 
here for fermions will also work in the bosonic case. 
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FIG. 4: (Color online) Average double occupancy dx/x during 
the expansion from a harmonic trap (Vtrap ~ O.OSt, C/init = 2t, 
A'' = 14, X = 2,3,...,7) with (a) a quench to U = 40t and 
(b) no quench. Dotted line in (a): dx^4{T — 0). Inset in (b): 
density profile at rt = 0, 15. (c) Entanglement entropy Sx 
for a: = 2, 3, 4 {U = 40t); horizontal lines are Sx{t = 0). 

tiferromagnetic MI with one particle per site. The idea 
would be to load such a BI in a trap for which that state 
is close to the ground state. This step is important in 
ensuring a low-entropy initial state that is in equilibrium 
[27]. One can then adiabatically reduce the strength of 
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ELECTRONIC PHYSICS AUXILIARY 
PUBLICATION SERVICE FOR: QUANTUM 
DISTILLATION: DYNAMICAL GENERATION 
OF LOW-ENTROPY STATES OF STRONGLY 
CORRELATED FERMIONS IN AN OPTICAL 
LATTICE 

Expansion from a box trap: dependence on the 
initial filling ninit 

In the main text, we have mostly discussed and pre- 
sented results for an initial filling of n-mit = 1-8 with 
N = 18. Here we present additional material for n-mit — 
1.2,1.4,1.6 (and, for comparison, the njnit = 1.8 data 
from the main text) to elucidate the dependence of the 
quantum distillation on Tiinit- Such data are presented 
for U/t = 8, 20, 40, 100 in Figs. 5-8, respectively. 

We find that for all ninit > 1 and all values of U con- 
sidered here, quantum distillation takes place. We distin- 
guish between the quantum distillation in a weak sense, 
namely for some sites (niir)) > {niir = 0)). Such be- 
havior is seen in the case of [/ = 8t (Fig. 5). The more 
pronounced effect of the formation of approximate Fock 
states, i.e., {ni{T)) w 2, is observed for large riinit and 
large U (see Figs. 6-8). 

Expansion from a box trap: dependence on the 
particle number at a fixed initial filling ninit 

We have also studied the dependence on the number 
of particles iV at a fixed riinit. The respective results are 
presented in Fig. 9 for U ~ lOOt and ninit = 1.8. Qual- 
itatively, the spatial extension and stability of the ap- 
proximate Fock states grows with N. The figure further 
shows [see Fig. 9(d)] that by plotting the data vs t/N, 
dx/x, calculated for different N, all results collapse onto 
essentially the same curve. This illustrates that the time 
for the formation of the quasi-Fock state scales linearly 
with the number of particles (or the number of doublons 
in the initial state, respectively). Most importantly, the 
qualitative behavior of the quantum distillation process 
is independent of A'^ at a fixed njnit 

For riiiiit < 1, no quantum distillation takes place at 
any U, and the dynamics in that case was studied in Ref. 
[8]. 



of U and an initial density of n = 1.8. In this quantity, 
the effect is visible for U >8t. 



Expansion from the harmonic trap: Reduction of 
the entanglement entropy 

Here we provide results for the reduction of the en- 
tanglement entropy for the expansion from a harmonic 
trap and the parameters of Fig. 4 of the main text 
(Vtrap O.OSt, N = 14, f/init = 2i), yet for two val- 
ues of the interaction U = AOt and U = IQQt after the 
quench. We have also performed runs with different dis- 
carded weights dp to check the quality of our results. The 
data presented in Fig. 11(b) suggest that the reduction of 
the entanglement entropy is (for x — 2,3 and U = lOOt) 

5,=2(t) < 0.15,=2(t = 0) (3) 
Sx^3{t) < 0.215,=3(t = 0). 

Note that we have not yet reached the minimum in the 
time-dependent reduction of Sx{t). At long times, the 
simulation is slowed down due to the relatively large en- 
tanglement of blocks with x >7. The figure unveils that 
the more accurate the simulation is {i.e., the smaller the 
discarded weight), the better the efhciency of the quan- 
tum distillation in reducing Sx{t) is captured. 

The last figure. Fig. 12, illustrates the dependence of 
the quantum distillation on first, U [panel (a)] and sec- 
ond, the effective density p = N^yVu-ap [panel (b)]. We 
use the parameters of Fig. 4 from the main text, except 
for varying U in Fig. 12(a) and varying TV in Fig. 12(b). 
We plot the average double occupancy and we observe 
that similar to the case of a box trap that was discussed 
in the main text, quantum distillation works as long as 
some doublons (signified by sites with (71^) > 1 at r = O) 
are present in the initial state. Also, as expected, the fi- 
nal U needs to be larger than 4i to induce the formation 
of a metastable state with doublons grouping together in 
the leftmost sites. 



Radius of the double occupancies 

The quantum distillation effect can be nicely illus- 
trated by plotting the radius Rd of the double occupan- 
cies: 

(2) 

where A^^ = X^i^i (^i) the total number of double occu- 
pancies. This radius is shown in Fig. 10 for several values 
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FIG. 5: (Color online) Average double occupancy dx{T)/x 
for a: = 2,..., 10 (d^r) = Er=o('^i(^)» and U = 8t: (a) 
riinit = 1-2; (b) ninit = 1.4 (c) ninit = 1.6; (d) ninit = 1.8. 
(time step O.Ol/t, discarded weight 10"**). 



FIG. 7: (Color online) Average double occupancy dx{T)/x 
for X- = 2,...,10 (d^r) = ELo(^~'W)) and U = 40f: (a) 

Tlinit = 1.2; (b) Tlinit = 1.4 (c) Tlinit = 1.6; (d) riinit = 1.8. 

(time step 0.01/t, discarded weight 10"'). 
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FIG. 6: (Color online) Average double occupancy dx{T)/x 
for x = 2,..., 10 (dxir) = ELo(^^~'W» and U = 20t: (a) 
ninit = 1.2; (b) riinit = 1.4 (c) riinit = 1.6; (d) ninit = 1.8. 
(time step 0.01//;, discarded weight 10~*). 



FIG. 8: (Color online) Average double occupancy dx{T)/x 
for a; = 2, . . . , 10 (d,(r) = E"=oK('^))) and [/ = 100/: (a) 
riinit = 1.2; (b) riinit = 1.4 (c) riinit = 1.6; (d) riinit = 1.8. 
(time step 0.01//, discarded weight 10~'). 



7 



X 

T3 



X 

•a 




10 20 30 
time 1 1 



time T t/N 



FIG. 9: (Color online) Dependence of the average double oc- 
cupancy dx{T)/x on the number of particles at a fixed filling 
of ninit = 1.8 {U = loot): (a) N = 18, L = 40 {Srt = 0.01, 
5p = 10"''); (b) TV = 36, L = 60 {Srt = 0.01, 5p = W'^); 
(c) iV = 36, i = 60 {Srt = 0.02, 5p = lO"'^); (d) d2(r)/2 
vs. time r/A'' (i.e., rescaled on the number of particles). 
The same qualitative behavior is observed, independently of 
N. The time for building up the quasi-Fock state is simply 
proportional to the total number of particles. Slight difi'er- 
ences between the results for different A'' at short times are 
due to differences in the initial state (i.e., boundary and A'^- 
dependent oscillations in the density profile). 
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FIG. 11: (Color online) Reduction of the entanglement en- 
tropy Sx{t) for X = 2,3,4,7 (bottom to top) and the pa- 
rameters of Fig. 4 of the main text (Krap = 0.05t, A'' — 14, 
;7init=2t): (a) U = 40i; (b) U = WOt. The runs were per- 
formed with a time step of 0.01/t and different discarded 
weights 5p (see the figure's legend). 




FIG. 12: (Color online) Average double occupancy ^2/2 for 
the set-up of Fig. 4 from the main text (Krap = O.OSt, [/init = 
2t) and (a) U/t = 2,8,20,40,100 (A^ = 14) and (b) A' = 
10, 12,14 ({/ = 40t). 



